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We study a static, spherically symmetric and asymptotic flat spacetime, assuming the hypersurface 
orthogonal Einstein-aether theory with an ultraviolet modification motivated by the Hofava-Lifshitz 
theory, which is composed of the z = 2 Lifshitz scaling terms such as scalar combinations of a 
three-Ricci curvature and the acceleration of the aether field. For the case with the quartic term 
of the acceleration of the aether field, we obtain a two-parameter family of black hole solutions, 
which possess a regular universal horizon. While, if three-Ricci curvature squared term is joined in 
ultraviolet modification, we find a solution with a thunderbolt singularity such that the universal 
horizon turns to be a spacelike singularity. 

PACS numbers: 04.25.dg, 04.50.Kd, 04.70.Dy 


I. INTRODUCTION 

Spacetime singularity is unavoidable in general 
relativity[l], which means the breakdown of our stan¬ 
dard theory of gravity in ultraviolet region. To establish 
the fundamental gravitational theory beyond general rel¬ 
ativity is one of the most intriguing question of physics. 
One may expect that this difficulty can be resolved by 
considering the quantum effect of gravity. However, un¬ 
fortunately, the perturbative quantization approach of 
general relativity losses the renormalizability unlike the 
other fundamental interactions. In other words, there ap¬ 
pears infinite numbers of divergent Feynman diagrams, 
and thus, the infinite counter terms are required to reg¬ 
ularize the gravitational quantum effects. Hence one 
way to quantize gravity is non-perturbative approach 
such as the loop quantum gravity [2] or the dynamical 
triangulation[3]. 

Another way is to find a new renormalizable gravi¬ 
tational theory. String theory [4] can be such a candi¬ 
date, but it has not been completed. Recently, Hofava 
proposed a gravitational theory with Lifshitz scaling [5], 
which is anisotropic scaling between space and time, i.e., 
t —>■ b~^t, x'’ —>■ b~^x''. This scaling defines a scaling 
dimension [t] = —z with [x*] = —1, which is restored to 
the ordinary mass dimension when z = 1. It is found 
that if we set z to the number of the spatial dimension, 
the dimension of the gravitational constant becomes zero, 
which means the gravitational force acquires renormaliz- 
ability at least at a power-counting level. This gravita¬ 
tional theory is called Hofava-Lifshitz (HL) theory whose 
action includes higher spatial curvature terms up to cubic 
order as counter terms of renormalization. 

Although there is a need for further investigation to 
confirm whether HL theory is truly renormalizable or 
not [6], strong gravitational phenomena, such as cos¬ 
mological singularity avoidance [7] and the black hole 
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solution[8] have been studied by several authors. In par¬ 
ticular, study regarding the spacetime structure is par¬ 
ticularly intriguing frontier. Since the spacetime in HL 
theory losses local Lorentz symmetry due to Lifshitz scal¬ 
ing, the causal structure is drastically changed. If there 
is Lifshitz scaling with z ^ 1, the dispersion relation of 
the signal particle is modified as ~ and then 
the sound speed is given by c ^ k^~^. In consequence, 
the sound speed almost diverges if the particle is in an 
extremely high energetic state. One may consider it is 
impossible to define the casual horizon due to such an 
instantaneously propagating particle. However it is not 
the case. In the context of the Einstein-aether (®-) the¬ 
ory [9] which has some equivalence to the infrared limit of 
HL theory [10], there still exists a causal horizon for such 
an extreme energetic particle. Although ae-theory itself 
is a toy model as a Lorentz violating gravitational theory, 
it is found that the HL theory is reduced to as-theory in 
infrared limit if the aether is restricted to be hypersurface 
orthogonal. In other words, the aether Up, is constrained 
to a gradient of some scalar field cp a.s oc 

If we set ip to time variable t, the HL action in ADM 
formalism is restored. Therefore p is called a “khronon. 
In this theory, the spacetime is expressed by a series of 
three dimensional spacelike hypersurface p = constant, 
denoted by Particles with infinite propagating sound 
speed travel along to Then, if there is a spacetime 
structure such that E^^^ parallels to a timelike Killing vec¬ 
tor namely u-^ = 0, any particles cannot escape from 
this surface at least in a spherically symmetric spacetime. 
Therefore this surface is a static limit for such energetic 
particles and it is called universal horizon[ll, 12]. The 
properties of the universal horizon have been so far inves¬ 
tigated on the following subjects : a static and spherically 
symmetric exact solution in ^-theory with the universal 
horizon[13, 14], the existence of the universal horizon[15- 
17], a charged black hole solution [18, 19], thermody¬ 
namical aspects[13, 20, 21], a ray trajectory in a black 
hole spacetime [22], and a formation via gravitational 
collapse [23, 24]. 

Then, one might ask a question whether the univer¬ 
sal horizon exist even if the Lifshitz scaling terms such 
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as higher spatial curvatures are present. In this paper, 
we shall consider the backreaction to the black hole solu¬ 
tion in ae-theory by the Lifshitz scaling terms. In other 
words, we investigate the black hole solution and prop¬ 
erties of the universal horizon with the Lifshitz scaling 
terms, which is the ultraviolet modification of gravity. 
As a first step to clarify the effect of gravitational Lif¬ 
shitz scaling, we shall consider only the z = 2 scaling 
terms. 

This paper is organized as follows: In §11, the action we 
considered is shown. We include the z = 2 Lifshitz scal¬ 
ing terms such as the quadratic spatial curvature to the 
action of ae-theory with hypersurface orthogonal aether 
and give the basic equations. The propagating degree of 
freedoms in this theory, i.e., the graviton and the scalar- 
graviton are also discussed. After giving the set up of our 
static, spherically symmetric system in §III, we classify 
numerical solutions depending on the coupling constants 
of the Lifshitz scaling terms in §IV. We then discuss the 
properties of a black hole solution and a thunderbolt sin¬ 
gularity in §V. §VI is devoted to conclusion of this paper. 


II. iE-THEORY WITH LIFSHITZ SCALING 

A. non-projectable HL gravity v.s. the ae-theory 
with Lifshitz scaling 

In order to see the behavior of the black hole solution 
with the backreaction from Lifshitz scaling for 2^1, we 
shall consider non-projectable HL gravity theory. In [25], 
its most general action is given as 


aether field is described by a gradient of a khronon 
field (fi as 




(2.4) 


Since the lapse function N in the HL gravity theory is 
related to this khronon field (p in the ae-theory as N = 
exp[(^], the aether field corresponds to the gradient of a 
lapse function <I>* in the low energy IR limit. 

Hence we will study the ae-theory with additional Lif¬ 
shitz scaling terms in order to discuss black hole solutions 
in the non-projectable HL gravity. Although the singu¬ 
lar behavior on the Killing horizon can be avoided by 
adopting Painleve-Gullstrand coordinate, it must be sin¬ 
gular at the universal horizon where the aether becomes 
normal to the timelike Killing vector. On the universal 
horizon, the khronon field p diverges, in other words, 
there is no continuous time coordinate beyond this hori¬ 
zon in the (3-|-l)-decomposition. We then reformulate 
the theory in covariant manner rather than the ADM 
(3-|-I)-decomposition. It is of great use to avoid coordi¬ 
nate singularities. 

In this paper, as a first step, we shall restrict the ul¬ 
traviolet modification terms only to simple scalar terms 
with z = 2 scaling such as TZ^ instead of considering all 
possible terms. The reason why only scalar terms are 
included is that those terms give the dependence in 
the dispersion relation of the scalar-graviton. Then one 
expects that the property of the horizon for the scalar- 
graviton which is generally singular in m-theory will be 
drastically altered. It would be appropriate terms to see 
the backreaction effect by the Lifshitz scaling. 


Ihl = J dt(f xN^/^ {Ck + Cp) (2.1) 

with 

Ck := a - A/C^) 

Cp := - (V.=i + V.=2 + V,=3) , (2.2) 

where N, g^^ij and K-ij are a lapse function, 3-metric, and 
an extrinsic curvature, respectively, and the potentials 
are defined by 

Vz=i := 7o'^ + 7i^*^' 

V.=2 := 73($*$*)" + •••+76(^^$*)7^-b•••+7lo7^^ 
V.=3 := 7ll(^^^*)^ + •••+7367^*^2?*P,7^. (2.3) 

a, A and 7 „ (n = 0,1,3--- ,36) are the coupling con¬ 
stants. The potential terms include not only the higher- 
order terms of the spatial curvatures TZij, TZ := TZ\ but 
also the non-linear terms of the gradient of a lapse func¬ 
tion <l)i := VilnN. Vz =2 and V 2=3 consist of 8 and 26 
independent terms, respectively [25]. 

The IR limit of non-projectable HL gravity theory is 
equivalent to the ae-theory with a hypersurface orthogo¬ 
nality condition[10], from which a spacetime is foliated by 
three dimensional spacelike hypersurface Then the 


B. The action and disformal transformation 

We consider the Einstein-aether gravity theory with 
z = 2 Lifshitz scaling terms, which action is given by 

^ ^ + '^(uv)] 

£(ir) = , 

■^(uv) = + yu^TZ + g 2 TZ^) ■ (2-5) 

where i? is a four dimensional Ricci scalar curvature, G 
is a gravitational constant mpi := Xj'/G is a Planck 
mass, which may corresponds to a typical Lorentz vio¬ 
lating scale, and the aether field is a dynamical unit 
timelike vector field. is defined by 

MZp ■■= Cl3^'>^"a + C25f^J''p - ClAuVg^p , 

( 2 . 6 ) 


^ Note that the Newton gravitational constant Gn is different from 
the gravitational constant G appeared in the action. Taking the 
weak field limit, we find their relation (3.11). 
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where C13 := Ci + C3, C14 := Ci + C4 with {ci} {i = 1-4) 
being the coupling constants in the ^-theory. 

The ^-theory given by /1 (ir) with a hypersurface or¬ 
thogonality condition is equivalent to the IR limit of non- 
projectable HL gravity theory [10]. The relation between 
both coupling constants is given by 

1 - Ci3 1 - C 2 

To 1 

a 


1 - Ci3 ’ 

7l Ci4 

a 


1 ’ 1 ■ (2.7) 

1 - Ci3 a 1 - Ci3 

A covariantized three dimensional Ricci scalar curva¬ 
ture TZ is given by 

7 ^ = R-{V^up){V^u°^) + {Vc,u°‘f 

-f 2Va[u“-u“(V;3M'^)], (2.8) 

where := is an acceleration of the aether 

corresponding to and iP' := . 

£(uv) is introduced as an ultraviolet modification mo¬ 
tivated by HL theory with z = 2, which is composed of 
the scalar combination of TZ and iP. The coupling con¬ 
stants in the z = 2 Lifshitz scaling are rewritten as 

I3i := 167173 , (32 ■■= 167476 , 32 := 1674710 . (2.9) 

We ignore the other coupling constants, i.e., 74 = 75 = 
76 = 77 = 78 = 79 = 0. 

Performing quadratic order perturbation of the action 
(2.5) around Minkowski spacetime, we find two types 
of the propagating degree of freedom. One is a usual 
helicity-2 polarization which corresponds to the gravi¬ 
ton. The other is helicity-0 polarization what we shall 
refer to as a scalar-graviton. The dispersion relations are 
given by 


Wg = 


wl = 


1 


1 - Ci3 


( 2 . 10 ) 


(Cl3 + C2)(2 — C14) ^^2 

Ci4(l — Ci3)(2 -|- Ci3 -I- 3c2) 


This transformation (2.14) means a rescaling of timelike 
separation between two spacelike hypersurfaces with fix¬ 
ing three-dimensional space. Under the transformation, 
the action is invariant if each coupling constant changes 
as 

Ci3 - 1 = ct(ci 3 - 1) , Ci3 -I- C2 = Cr(ci3 -|- C2) , 

Ci4 = C14 , 52 = 52 , /3l = /?! , ^2- (2.15) 

Remarkably, the coefficients proportional to k"^ in (2.10) 
and (2.11) are changed to times after the transfor¬ 
mation, whereas, terms are invariant. This means the 
propagating speeds of each gravitons in infrared limit are 
scaled as a~^P but that of the graviton in ultraviolet 
limit is unchanged. This property holds even if the all 
possible higher curvature terms motivated by HL theory 
are considered (see Appendix.B). 


C. The basic equations 

To derive the basic equations, we shall start from tak¬ 
ing the variation of action (2.5) with respect to and 
i.e.. 


(5/ = 


: J ■ Sg^^'' + 2JE^ ■ Su^^] (2.16) 


16nG , 

where, 

E,. := Ej^^^-m;,^[g2EliP+P^Eli^^+f32Elfp^], 

(2.17) 

iE^ := 4®) - m-^[g2MlsP + /3i7Ef) + ^2^^^^^], 

(2.18) 

The infrared portions, E^!^^ and are defined by ^ 

E^i, ^ .— G~J yg (VaU^^ 9 


(—)'■ 

( 2 . 11 ) 

+Vq, 


\m^ij 





Note that the infrared portions which is proportional to 
correspond to ae-theory’s one[26]. 

As is the case of the ®-theory[27], we find an invari¬ 
ance in the above model under the following disformal 
transformation; 

g^iu = gfj,i> + (1 - a-)u^u^ , u'^ = , ( 2 . 12 ) 

where cr > 0. This transformation can be simplified by 
introducing three metric on the spacetime hypersurface 
which is defined by 

7m^ 9iiv E u^Ui /. (2.13) 

The disformal transformation (2.12) is rewritten as 

liiu = Itiu , . (2.14) 


(2.19) 

:= + ci 4 ha (V^u“) , ( 2 . 20 ) 

where G^y := — ^Rg^i, is Einstein tensor and 

:= (V«u'5) 

= Ci3V,u^ + C2(V„u“)<5't - cuu>^u, .(2.21) 

The ultraviolet portions, Elfp\ Ejfp\ 

and are obtained 


as 


^ The round and square brackets in the tensoral index are a sym- 
metrization and anti-symmetrization symbols, respectively, i.e., 


(2.14) 
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= TZRfj,^ - + {"^‘^Tl)g^.u + '^a{a°‘TZ)g^^ - V(^V^,)7^ 

+Va [%-«l,(V“7^) - M“(V(^-«^.))7^ - ■U(^(VJ,)u“)7^ + (V“M(^)u^)7^ - 2-«“■U(^(VJ,)7^)] , (2.22) 


£’(/3i) 

r(/32) 


^032) 

H' 


2 ^ 4” ^a [^416 u ‘2ii u 

1 


--ii^Tlg^,^ - nUfj,Uu + u^Tlf,^ - V(^Vj,)M^ + V. 


-u^u“(V • + Wp{u'^u°‘u^)gfj,^ + u^M“V(^Uy) + u^M(^Vy)WQ 


(2.23) 

nu°‘Uf,u^ + V“(m^ 7 ^i,) 

- m^(V“m(^)mj.) - 2 V(^(mi.)U^m“) , 


(2.24) 

(VaV^w“ - V/.Vcw“)7^ + (V«u“)(V;.7^) - (V;.M“)(Va7^), (2.25) 

2 m2m^(V^u“)-2V„(m2u^u“), (2.26) 

7^Ma(V^u“) - V;3(MM«^7^) + u2(V„V^m“) - V„(u2V^m“) + V„V^(w2m“) + V^(w2v • m) 
-u2V^(V-u)-V^V„(m2u“). (2.27) 


Note that = 0 and = 0 are not the basic equa¬ 
tions, because the constraint of the aether field has not 
been taken into account. To find the basic equations, we 
usually have to introduce a Lagrange multiplier. Instead 
expressing the aether by /\/—UaU°‘ where is 

an arbitrary timelike vector field, we find the basic equa¬ 
tions from a variation of U^. Since a variation of is 
given by 


6u^^ 


--U>^UaU0{Sg“^) 


-f 


(( 5 ^ 0 , -I- M^Mg) 




we find the basic equations as 


(TEf^ii '^Ufj^Ui/ — 0 , 

(f?/ia 4 “ U^Uq,) 7 E — 0 , 

MqU“ = — 1. 


(2.28) 


(2.29) 

(2.30) 

(2.31) 


Here we rewrite the basic equations in terms of the aether 
field with the normalization condition (2.31). 

If the aether field is hypersurface orthogonal as we have 
assumed here, we can take a variation with respect to the 
khronon field ip instead of U^. Since the aether field 
is given by Eq. (2.4), the variation of p is found by using 
the relation: 


and 


^-(V°^p)(Vcp) 


(2.34) 


with the definition (2.4). 

In the case of the hypersurface orthogonal aether field, 
although the basic equations are given by Eqs. (2.33) 
and (2.34) with (2.4), these equations contain higher- 
derivative terms of the khronon field p. If spacetime is 
static and spherically symmetric, however, we may find 
simpler equations, which are the original basic equations 
(2.29)-(2.31). It is because the hypersurface orthogo¬ 
nality of the aether is automatically satisfied for spher¬ 
ically symmetric spacetime, and then the original basic 
equations are reduced to the basic equations with hyper¬ 
surface orthogonality Although those equations are 
equivalent, Eqs. (2.29)-(2.31) are written in terms of the 
aether field u^, then those are the second-order differ¬ 
ential equations of u^. For this reason, we shall adopt 
the (2.29)-(2.3I) as the basic equations in the rest of this 
paper. 


III. SPHERICALLY SYMMETRIC “BLACK 
HOLE”: SET UP 


Su^ = 


/3) 

( 5 '"“ +«''«“) 


yJ-g°‘^{Vc,p){yi3p) 


(Sg^^) 

V„(M. (2.32) 


We discuss a static and spherically symmetric space- 
time with asymptotically flatness. In order to avoid a 
coordinate singularity at horizon, we adopt the following 


The resultant basic equations are 

T u^Ui/(^J^Q^u ) “h = 0 , 


^ The equality of these set of equations holds if the spacetime is 
regular everywhere. This can be proven by considering volume 
(2.33) integral of (2.34) and using Gauss’s theorem[28]. 
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metric ansatz like the Eddington-Finkelstein type: 

ds^ = —T{r)dv^ + 2B{r)dvdr + r^dfl^ , (3.1) 

where v is an ingoing null coordinate and B > 0. The 
aether field in this coordinate system is assumed to be 

u>^ = (a{r),b{r),0,0j , (3.2) 

where the function 6 (r) is fixed by the normalization con¬ 
dition (2.31) as 


^ a(r)^T(r) - 1 
^ ^ 2a{r)B{r) 


(3.3) 


In this spacetime, there exists a timelike Killing vector 
:= ( 1 , 0 , 0 , 0 ) associated with the time translational 
invariance. 

Since the basic equations (2.29)-(2.31) in this ansatz 
take quite complicated form, we omit to show it explic¬ 
itly. Instead, the structure of the basic equation is illus¬ 
trated. Substituting (3.1) into the basic equation (2.29)- 
(2.30), we find there are five non-trivial and independent 
set of equations : {v,v), (f,r), (r,r) and (0,0) compo¬ 
nents of (2.29) and projection of (2.30), where is a 
“radial” spacelike unit vector perpendicular to u^. From 
the discussion in [29], we find following two constraint 
equations : 


C^ = 0, C"' = 0, (3.4) 

where is defined by 

C^^ _ 2 {JE^u°‘) . (3.5) 


These equations include one fewer r derivatives than the 
rest portion of the basic equations, and they are auto¬ 
matically preserved by solving the other equations with 
respect to r-evolution if (3.4) are satisfied on an “initial” 
constant r-surface. The rest of the equations, namely 
{v,v) and {6,9) components of (2.29) and component 
of (2.30) give the evolution equations with respect to 
T{r), B{r) and a{r). 


A. Asymptotic behaviour 


Since we assume an asymptotic flatness, the asymp¬ 
totic values of the variables are given by 

T{r) 1, B{r) 1, a(r) ^ 1, h{r) 0. (3.6) 


In order to investigate the asymptotic behavior of the 
solution, we perform an asymptotic expansion around 
Minkowski spacetime, that is, the functions T(r), B{r) 
and a{r) are expanded as a series of l/r as 


T{r) 

B{r) 


Ti 

= 1 + — 
r 




Bi B 2 

= 1 + —+ — 


r 

ai 


a{r) = 1 + —+ — + — + —+ 


a2 

2 


^3 ~ 

El 

,3 


r 

03 

.3 


r4 

~ 4 

a4, 

7^ 


(3.7) 


Substituting these series into the basic equations, and 
solving them order by order, we find the expansion coef¬ 
ficients as 


J 


Ti 


T 4 


arbitrary, Tz = 0 , Tg = , 

{(4ci4 + 19)ci4 - 54ci 3}T4 + 192m-i2(ci452 - P2)T^ + 48(ci4 - 2ci3)(4a2 - 3T2)a2 

192(2-C 14 ) 


Bi = 0 , B 2 = 


ciiT^ 

16 


B. = - 


C 14 T 1 


12 


B. = 


3ci4{cli + 14ci4 - 36ci 3 + 4)T4 + 256m-i2(2ci4 - l)(ci452 - P2)T^ + 192(1 - 013 ) 01402(402 - ST^) 

512(ci4-2) 


oi = —^ , 02 = arbitrary, 03 = — 


^14 ^rp3 I rp 

-^T,+T,a2 


1 


CI 4 = 


1920(ci4 — 2)ci23 


[5c2{5ci4(2ci 4 — 1) 24} H- 18ci4(ci4 — 2) H- ci3{32c^4 + llci4 H- 30(4 — 9c2)} — 270cf^]Ti 


+ 48[[(2 — Ci4)ci4 + 10c2(ci4 — 5) + Ci3{30Ci 23 + Ci4(ci4 + 8) — 50}]a2 + 20ci23^pi^(ci452 — P2)]Ti 


1 + Ci3 — Ci4 

2 — Ci4 


a 


2 

2 • 














6 


The asymptotic behavior of the function b{r) is given by 


1. horizons in the infrared limit 


with 

3Tf 

02 = 02 -^ , 

ci4(3ci23 + 2c 2 + 2)(3Tf — 8 a 2 )T'j 

-. (3.9) 

The important point is every order of these functions, 
at least up to the eighth order, is described only by two 
arbitrary coefficients, Ti and 02 as in the case of the 
ae-theory[30]. Additionally, the effect of the 32 and /32 
terms, that is, the contribution from the fourth spatial 
derivative terms first appears in the fourth order coef¬ 
ficients and 04 . /3i appears after the fifth order 

coefficients, which we have not shown here because they 
are so lengthy. 

The free parameter Ti is accosted with a black holes 
mass. From the discussion in [31-33], the black hole mass 
M as a Noether charge with respect to time translational 
symmetry is given by 


M = - 


II 

2G 



Ti 

2Gn ’ 


(3.10) 


where 

Gn:=g(i-^)'' (3.11) 

is the observed Newton constant. 

The parameter 02 can be fixed from the analyticity of a 
black hole horizon for the scalar-graviton in the infrared 
limit, but it becomes a free parameter when we include 
the z = 2 Lifshitz scaling terms as will be discussed later. 
This free parameter 02 may characterize the distribution 
of an aether cloud around a black hole. 


B. Black hole horizons 

In the ae-theory or the HL gravity theory, the metric 
horizon, which is the r-constant null surface of (3.1), is 
not generally an event horizon. As shown in Eqs. (2.10) 
and ( 2 . 11 ), the sound speeds of the graviton and the 
scalar-graviton depend on the coupling constants. As 
a result, without tuning of the couplings, they generally 
differ from unity. The metric horizon only means a static 
limit for a propagating mode with the sound speed being 
unity. For this reason, we first reconsider the horizons of 
the aether black hole. 


4 We shall refer Mpi := I/^/Gn as a observed Planck mass which 
is related to the observed Newton constant Gn- Note that the 
Planck mass mpi := 1/\/G appeared in (2.5) is rather than re¬ 
lated to the Lorentz violating scale. 


Firstly, we shall consider the low-energy infrared limit 
of the graviton and scalar-graviton. Since the relevant 
parts in ( 2 . 10 ) and ( 2 . 11 ) are terms, the sound speed 
of the graviton cg and that of the scalar-graviton cs in 
the infrared limit are given by 


Cg 


1 — C13 ’ 


Cs 


(Cl 3 + C2)(2 — C14) 
Cl 4 (l — C 13)(2 -I- C13 -I- 3c2) 


.(3.12) 


Note that under the transformation (2.14), each sound 
speed is changed as 


Cg = O' ^Cq, and Cg = a ^Cg . (3.13) 


Therefore, if we set ct = Cq or Cg, we find a frame in 

which either the sound speed of the graviton or that of 

the scalar-graviton is unity. Thus, we can adjust the 
horizon for the graviton or that of the scalar-graviton to 
the metric horizon by an appropriate disformal transfor¬ 
mation. Explicitly, by performing the foilwing disformal 
transformations; 

9\?l ■= + (1 - CG)Uf,u^ , (3.14) 

9 III ■= 9n’' + (1 - 4)unUo ■ (3.15) 


the graviton horizon or the scalar-graviton horizon are 
located on the r-constant null surfaces of the effective 
metrics (3.14) and (3.15), respectively. In the Eddington- 
Finkelstein ansatz (3.1), the r-constant null surfaces of 
the graviton and the scalar-graviton are given by 


TG(r)=0, rs(r)=0. 


respectively, where 


(3.16) 


Tg ■■= 

Ts ■■= -9ii 


T- (1 
T- (1 


1 + a^T 


4 ) 


2a 

1 + a^T 
2a 


(3.17) 

(3.18) 


respectively. In Appendix.A, we present the transforma¬ 
tion of Eddington-Finkelstein type metric (3.1) under the 
disformal transformation (2.14). 


2. horizons in the ultraviolet region 


In turn, we shall focus on the propagation of the gravi¬ 
ton and the scalar-graviton in the high energy limit. Al¬ 
though the sound speed of the graviton is the same as 
that in the infrared limit, the sound speed of the scalar- 
graviton in the high energy limit turns to be 


c 


2 

s 


8 (ci3 -I- C2)g2 / k \ 
2 -I- Ci3 -I- 3 c2 \mpi) 


(3.19) 
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which depends on the three momentum k. Thus, the 
sound speed can increase to infinitely high in an ulti¬ 
mately excited state. In this situation, the r-constant 
null surface given by Eq. (3.15) is no longer an event 
horizon. An event horizon of a black hole must be the sur¬ 
face whose outside region is causally disconnected from 
the inside for any propagation modes even with an in¬ 
finite sound speed. Otherwise, an inside singularity is 
exposed. 

The above case can be resolved to consider the spe¬ 
cial aether configuration. The ultimately excited scalar- 
graviton should propagate along the three dimensional 
spacelike hypersurface. In other words, any future di¬ 
rected signal must not propagate against the direction 
which Lp decreases. Thus, such an excitation mode must 
be trapped inside a surface where the hypersurface 
is parallel to the timelike Killing vector namely, u ■ ^ 
vanishes. This is the concept of the universal horizon 
which is regarded as a real black hole horizon in Lorentz 
violating spacetime[II]. 


IV. SPHERICALLY SYMMETRIC “BLACK 
HOLE”: SOLUTIONS 

To find a black hole solution with the z = 2 Lifshitz 
scaling terms, we shall solve the basic equations numer¬ 
ically. Our strategy is as follows, (i) To impose the 
boundary conditions near the asymptotically flat region 
by applying (3.8), that is, to give “initial” values of the 
variables T{r),B[r), and a{r) and their derivatives at 
infinity, (ii) To integrate from an appropriate distant 
spatial point toward the center of a spherical object. 


A. black hole solution in the infrared limit : 

the case of 32 = di = /32 = 0 

First we show the result for the case of the ae-theory, 
i.e., g 2 = I3i = P 2 = 0. It gives a black hole solution 
in the low-energy infrared limit. The numerical black 
hole solution is shown in FIG.I, which was already found 
in [30], for the coupling constants C 13 = 0.100,C 2 = 
—6.135 X 10 “"^,ci 4 = 0.100. In this solution, the black 
hole mass is chosen as GnM = 0.5. Note that our unit 
is fixed by setting Ti = —1, so that the normalization 
length is tm := 2 G^M = 1 . 

Since, the scalar-graviton sound speed is set to unity, 
the scalar-graviton horizon coincides with the metric 
horizon : T(rsH) = 0. We find rsn = 1.010 = 2.02GnAI, 
which is a little larger than the Schwarzschild radius. 

The graviton horizon locates inside that of the scalar- 
graviton, i.e., rcH = 0.998 < rsn- It is because the 
graviton sound speed is faster than the scalar-graviton’s 
one. The most outer universal horizon is formed inside 
these two horizons, i.e., run = 0.720 < t-ghi^'sh- Addi¬ 
tionally, more than one inner universal horizons is formed 
due to the rapid oscillatory behavior of the function U (r) 



FIG. 1: Black hole in the ae-theory[30]. The coupling con¬ 
stants are set to C13 = 0.100, C2 = —6.135 x 10~^, C14 = 0.100. 
We tune the parameter as 02 = 1.112 x 10“^, which gives a 
regular scalar-graviton horizon. The graviton sound speed 
and the scalar-graviton’s one are Cq — 1.111 and c| = 1.000. 
The black hole mass is normalized as GnM = 0.5. 

The blue, red and green curves indicate the functions 
T{r),B{r) and a(r), respectively. The scalar-graviton 
horizon is given by rsn = 1.010. The dashed purple curve 
indicates Tq, from which we find the graviton horizon as 
rcH = 0.998. The dashed brown curve indicates U{r) := u-G 
which zero point gives the position of the universal horizons. 
The outermost universal horizon radius is run = 0.720. 


near the central singularity. It means, this solution has 
a causally disconnected region for low energy particles 
even if the particle with z ^ \ Lifshitz scaling is taken 
into account In this sense, this solution is regarded as 
a black hole in the low-energy infrared limit. 

The important point of this solution is that a physi¬ 
cal singularity generally appears on the scalar-graviton 
horizon (= the metric horizon in the present case), if we 
do not tune the parameter 02 . In order to regularize 
the scalar-graviton horizon, we must choose an appropri¬ 
ate value for the parameter 02 as a boundary condition, 
which is 02 = 1.112 x 10“^. Then, the function 
which is the “coefficient” of l/Ts in the evolution equa¬ 
tion of B[r), must vanish on the scalar-graviton horizon 
T = ?'SH- We present the detailed analysis of the regu¬ 
larity on the horizons in Appendix C. As a result, for 
a regular ae-black hole solution, there remains only one 
free parameter Ti, which is associated with the black hole 
mass, just as is the case of the Schwarzschild solution in 
general relativity. 


® Note that an instantaneous propagating mode appears when the 
interaction between khronon and matter field is taken into ac¬ 
count even if the higher spatial derivative terms in action are 
absent[12, 34]. In our discussion, however, we focus only on the 
gravitational part of the theory whose action is given by (2.5) 
without T(uv). 











B. black hole solutions with Lifshitz scaling : 

the case of / 3 i 7^ 0 and 32 = /92 = 0 

When the higher-order aether correction is taken 
into account, i.e. j3i the solution turns to depend on 
02 as well as Ti unlike ae-black holes. Therefore one may 
consider the ultraviolet correction does cure the singu¬ 
lar behavior on the scalar-graviton horizon appeared for 
the infrared-limit theory. In fact, if we assume only 
term {g 2 = /?2 = 0 ), we find that there is no singular be¬ 
havior on any horizon in general. The detailed discussion 
is developed in Appendix. C. 

In this subsection, we consider the case of 7 ^ 0 with 
92 = 1^2 = 0. There are five types of “black hole” solu¬ 
tions, which are classified in two dimensional parameter 
space, namely, in ( 02 ,/3i) plane. Note that we use the 
unit of Ti = —1. We give a classification of these solu¬ 
tions in TABLE I, which phase diagram of these solutions 
is shown in Fig. 2. 


solution 

region 

horizons 

singularity 

BH 

iBH(i) 

I(i) 

rcH < rsH & no run 

r = 0 

A 

iBH (ii) 

I(ii) 

rcH < rsH & no run 

0 < r < rcH 

A 

uBH (i) 

II(i) 

ruH < rcH < rsH 

r = 0 

0 

uBH (ii) 

Il(ii) 

ruH < rcH < rsH 

0 < r < ruH 

0 

iTS 

III 

no horizons 

r = rsH 

X 

tNS 

IV 

no horizons 

r = 0 

X 


TABLE I: The classification of the solution for the case of 
/3i 7 ^ 0 with P 2 = 92 = 0. The region I-IV correspond to the 
areas shown in Fig. 2. rsH,rGH and run are the positions 
of the scalar-graviton horizon, the graviton horizon, and the 
universal horizon, respectively. A means that the solution 
describes a black hole for gravitons and scalar-gravitons, but 
may becomes naked for high-energetic Lifshitz scaling test 
particles with 2 > 1. 

We explain each solution in due order: 

(1) iJBH (i) [an infrared black hole with a central singu¬ 
larity] : If Pi and 02 are set to be in the light blue colored 
region in Fig. 2 (region I(i)), we find a kind of black hole 
which possesses the graviton horizons. The typical nu¬ 
merical example is shown in Fig. 3(a). There are gravi¬ 
ton and scalar-graviton horizons. In the present case, 
the dispersion relations of graviton and scalar-graviton 
are given by ~ which means these horizons co¬ 
incide with ae-theory’s ones. Hence the gravitons and 
scalar-gravitons cannot escape from the inside of these 
horizons. For low-energy particles, they also play a role 
of horizon too. In this sense, one may regard this solution 
as a type of black hole, which we call an infrared black 
hole {iBH). 

The spacetime and aether field are regular except at 
the center. However, since there exists no universal hori¬ 
zon, this solution has no causally disconnected region. 
The universal horizon turns to be genuine causal bound¬ 
ary due to the Lifshitz scaling with 27 ^!. Therefore the 
singularity at the center is exposed if non-gravitational 



(a) The phase diagram of the solutions 



(b) The enlarged phase diagram near the se-black hole 

FIG. 2: The phase diagram of the solution in the ( 02 ,/?i) pa¬ 
rameter plane, where The coupling constants are chosen as 
the same as Fig.l. The red dashed line which is /3i = 0 indi¬ 
cates the case of *-theory. The ae-black hole with a regular 
scalar-graviton horizon is shown by the red circle. The gen¬ 
uine black hole solutions which possess the universal horizon 
are discovered in the region II(i) and Il(ii). 


propagating modes with the 2 > 1 Lifshitz scaling are 
taken into account. Thus we conclude that this solution 
does not describes a true black hole in the strict sense but 
a type of naked singularity even if the graviton horizons 
exist. 

To clarify this situation, we shall depict the spacetime 
structure. Note that the Carter-Penrose diagram itself 
does not describe the causal structure of the solution 
due to the lack of Lorentz invariance. However, since it 
would be useful to understand the spacetime structure, 
we will show it for this solution. In Fig. 3(b), we 
illustrate the Carter-Penrose diagram for this solution, 
in which null rays propagate on ±45o direction. The 
metric (and scalar-graviton) horizon, which is one of the 
horizons in the Carter-Penrose diagram, is a horizon 
for the 2 = 1 Lorentz invariant particles or for the 
low-energetic infrared particles. For the 2 > 1 Lifshitz 
scaling high-energetic particles, however, it is no longer 
horizon, but the spacelike universal horizon will take its 
place. 

(2) iBH (ii) [an infrared black hole with a singular 
spherical shell ]: This solution can be found in the deep 
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(a) The evolution of the each components of the metric and 
aether. 


r = 0 : singularity i~^ 



(b) The Carter-Penrose diagram 

FIG. 3: The typical example of iBH (i) (region I(i)). In 
the top figure (a), We choose {a 2 , Pi/m^i) = ( — 1.000,1.000). 
The remaining coupling constants and the boundary condi¬ 
tions are set to the same values as those of Fig. 1. The 
blue, red and green curves indicate the functions T{r), B{r) 
and a{r), respectively. The dashed purple and dashed brown 
curves indicate Tg and U{r) := u ■ The graviton horizon 
is found at roH = 1.030, while the scalar-graviton horizon 
exists at rsH = 1.120. Since U{r) does not vanish, there is no 
universal horizon, namely, no causal boundary. In the bot¬ 
tom figure (b), the dashed curves indicates r =constant sur¬ 
face, especially, the blue, purple and brown curves represent 
the scalar-graviton horizons, the graviton horizons and the 
universal horizons, respectively. The spacetime singularities 
are represented by the waving curves. The future (past) null 
infinity, the future (past) timelike infinity and the spacelike 
infinity are indicated by (J^~), (*~) and respec¬ 

tively. The ultimately excited particle with z > 1 Lifshitz 
scaling propagates along the ip =constant surface, which is 
indicated by the solid green curves. 


blue colored region in Fig.2 (region I(ii)). Although 
graviton and scalar-graviton horizons are formed, a 
singularity appears at r = rss > 0. As mentioned in (1), 
this singularity is not causally disconnected from infinity 
due to the absence of the universal horizon. Therefore, 
although this solution behaves as a black hole for gravi¬ 
tons and scalar-gravitons as well as Lorentz invariant 
z = 1 particles, it turns to be a naked singularity for 
high-energetic particles with the z > 1 Lifshitz scaling. 


In this sense, we also classify this solution as iBH. 

The difference from the case (1) is that the singularity 
shapes a spherical shell rather than a spacetime point 
with infinitesimal volume in iBH (i). In this paper, we 
shall refer it as a singular shell. 

To see the cause of this singularity, we shall focus on 
the structure of the evolution equation. We find that 
the evolution equation of B{r) which is a linear-order 
differential equation with respect to r (see Appendix C) 
turns to be singular at r = More specifically, the 
divergence of B'(r) results in this type of singularity. 
The typical numerical example and the Carter-Penrose 
diagram are shown in Fig. 4(a) and (b), respectively. 



(a) The evolution of the each components of the metric and 
aether. 



(b) The Carter-Penrose diagram 

FIG. 4: The typical example of iBH (ii) (region I(ii)), In 
the top figure(a). We choose, {a 2 , Pi/m^i) = (0,1.000). The 
remaining coupling constants and the boundary conditions 
are set to the same values as those of Fig. 1. The blue, red 
and green curves indicate the functions T{r),B{r) and a(r), 
respectively. The dashed purple and dashed brown curves 
indicate the Tq and U{r). The graviton horizon is found at 
rcH = 0.949, while the metric horizon exists at rsu = 0.970. 
The singularity appears at rss = 0.792, which gives the radius 
of the singular shell. In the bottom figure(b), the conformal 
structure is depicted. The meaning of the curves and symbols 
in this figure are same as those of FIG. 3(b). 

(3) uBH (i) [an ultimate black hole with a central 
singularity ]: The black hole solution which possesses the 
graviton, scalar-graviton and universal horizon and no 
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singularity except the center is found in the light orange 
colored region in Fig. 2 (region II(i)). Since the central 
singularity is hidden by the universal horizon, this type 
of the solution is a real black hole. Any particles with 
the z > 1 Lifshitz scaling as well as the Lorentz invariant 
z = 1 particles cannot escape from the inside of the 
universal horizon. We then call it an ultimate black hole 
{uBH). 

We shall show the typical example in Fig. 5(a). This 



(a) The evolution of the each components of the metric and 
aether. 


r = 0 : singularity 



(b) The Carter-Penrose diagram 

FIG. 5: The typical example of uBH (i) (region II(i)). In the 
top figure(a), we choose (o2,/3i/mpi) = (0.200,1.000). The 
remaining coupling constants and the boundary conditions 
are set to the same values as those of Fig. 1. The blue, red 
and green curves indicate the functions T{r), B{r) and a{r), 
respectively. The dashed purple and dashed brown curves 
indicate Tq and U{r). The graviton and the scalar-graviton 
horizons are found at tgh = 0.882 and rsu = 0.896, while 
the universal horizon exists at run = 0.320. Since there are 
a universal horizon and a central singularity, this solution is 
referred to a black hole solution. In the bottom figure(b), the 
conformal structure is depicted. The meaning of the curves 
and symbols in this figure are same as those of FIG. 3(b). 

solution is much similar to the ae-black hole illustrated 
in Fig. 1 except the oscillatory behavior of U{r) 
near the central singularity. Moreover, we also show 
the Carter-Penrose diagram of this solution in FIG. 5(b). 

(4) uBH (ii) [an ultimate black hole with a singu¬ 


lar spherical shell ]: When and /3i are in deep 
orange colored region in Fig.2 (region Il(ii)), we find 
the solution with graviton, scalar-graviton and universal 
horizon, however, there are a singular shell inside the 
universal horizon. Since the singular shell is covered 
by the universal horizon, any information from the 
singularity never be leaked into the outside. Thus we 
can regard this type of the solution as a real black hole. 
Then we also classify this solution as an ultimate black 
hole [uBH). 

The typical example and the Carter-Penrose diagram 
of this solution are shown in Fig. 6(a) and (b), respec¬ 
tively. 



aether. 



(b) The Carter-Penrose diagram 

FIG. 6: The typical example of uBH (ii) (region Il(ii)). In 
the top figure(a), we choose (o2,/Ii/”^pi) = (0.117,1.000). 
The remaining coupling constants and the boundary condi¬ 
tions are set to the same values as those of Fig. 1. The 
blue, red and green curves indicate the functions T{r), B{r) 
and a{r), respectively. The dashed purple and dashed brown 
curves indicate Tq and U{r). The graviton and scalar- 
graviton horizons are found at rcH = 0.972 and rsH = 0.983, 
while the universal horizon exists at run = 0.524. The sin¬ 
gularity appears at Vss = 0.392, which gives the radius of the 
singular shell. In the bottom figure(b), the conformal struc¬ 
ture is depicted. The meaning of the curves and symbols in 
this figure are same as those of FIG. 3(b). 


(5) iTS [an infrared thunderbold singularity]: For 
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the solutions in the light green colored region in Fig. 
2 (region III), a singularity always appears at the null 
Killing horizon of the scalar-graviton. On the singular 
shell, the function a(r), namely, the v component of the 
aether field diverges if we set 02 to be larger value than 
a critical one. Furthermore, it is found that the critical 
value of 02 (~ 0.377) which induces the a(r) divergence 
seems to be universal under any choice of j3i as shown 
in Fig. 2). We shall present one typical example and 
the corresponding Carter-Penrose diagram in Fig. 7. 
Near the spherical shell with the radius rsn = 0.865, 
the quadratic scalar of three-dimensional Ricci tensor 
TZfi^TZ^'^ diverges. Hence it is a physical singularity. 

One may wonder why such a singular behavior is oc¬ 
curred on the scalar-graviton horizon. More specifi¬ 
cally, from Appendix C, we have known that the ba¬ 
sic equations possess no dependence on the negative 
power of Ts unlike the ae-case. Therefore, the scalar- 
graviton horizon where Ts = 0 should be regular in 
general. To see in detail, we shall expand the evolu¬ 
tion equation of a(r) around rsn- From our numeri¬ 
cal analysis in Fig. 7, we find a“^(rsH -f e) ^ ^ and 
o^(^SH + e) S~^, where e <C 1 and (5^1. In Fig. 
7, we have check the divergence numerically at least for 
e > 1.000 X 10“^ and <5 >4.132 x 10“^. Near the singu¬ 
larity, where Ts(rsH -f e) := -C 1, the dominant term 
in the evolution equation of a{r) is given by 


a"{r) 


S(^Sh)2A,“ ■ 


(4.1) 


Thus, we find that a(r) shows singular behavior when 
A —>■ 0, where the limit gives the scalar-graviton horizon. 
Since this singularity is due to the aether field rather than 
the spacetime metric. 

This singularity seems to be closely analogous to the 
thunderbolt singularity. The thunderbolt singularity is 
proposed in the context of quantum black hole evapo¬ 
ration. In [35], the thunderbolt singularity is firstly in¬ 
vented in semi-classical analysis of (1 -I- l)-dimensional 
dilaton-coupled gravity with scalar field [36] which is 
renormalizable theory of quantum gravity. Furthermore, 
this type of singularity is also discovered in (1 -|- 1)- 
dimensional quantum field theory via complete quantized 
analysis [37]. 

In such a situation, a null singularity appears on the 
event horizon. Although a causally disconnected region is 
not formed because of the existence of a singularity, this 
singularity itself is not detected by any outside observer. 
As a result, it is not a naked singularity. They call it a 
thunderbolt singularity. 

In our case, the singularity of the aether field appears 
on a Killing horizon is null. Hence it behaves similar 
to a thunderbolt singularity for Lorentz invariant z = 1 
particles. We call it an infrared thunderbolt singularity 
(iTS). The thunderbolt composes of the singular aether 
field. 



(a) The evolution of each components of the metric and 
aether. 



(b) The corresponding evolutions of the curvatures. 



(c) The Carter-Penrose diagram 


FIG. 7; The typical example of iTS (region III). In the top 
figure(a) and the middle 6gure(b), we choose ( 02 ,/3i/mpi) = 
(0.500, 1.000). The remaining coupling constants and the 
boundary conditions are set to the same values as those of Fig. 
1. In the top figure (a), the blue, red and green curves indicate 
the functions T{r), B{r) and a{r), respectively. The dashed 
purple and dashed brown curves indicate the Tg and U{r). 
Near the spherical shell with the radius rsn = 0.865, the 
quadratic scalar of three-dimensional Ricci tensor 77^1.77,^'' 
diverges rather than four-dimensional Kretchmann invariant 
77^1,77^", which are denoted by the solid dark green curve and 
dot-dashed black curve in the middle figure(b), respectively. 
In the bottom figure(c), the conformal structure is depicted. 
The meaning of the curves and symbols in this figure are same 
as those of FIG. 3(b). 


(6) tNS [a timelike naked singularity without horizon]: 
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The singular shell without any horizon is found in the 
light red colored region in Fig. 2 (region IV), namely, 
/3i < 0. Since the singularity is timelike and there is no 
horizon, this type of the solution completely exposes its 
singularity. Note that that singularity is originated from 
the divergence of the evolution equation B{r) as is the 
case in iBH (ii) and uBH (ii). 

We shall show the typical example in Fig. 8 . 



(a) The evolution of each components of the metric and 
aether. 



(b) The Carter-Penrose diagram 


FIG. 8: The typical example of tNS (region IV). In the top 
figure (a), we choose ( 02 ,/3i/mpi) = (-1.000,-1.000). The 
remaining coupling constants and the boundary conditions 
are set to the same values as those of Fig. 1. The blue, red 
and green curves indicate the functions T{r),B{r) and a(r), 
respectively. The dashed purple and dashed brown curves 
indicate the Tq and U{r). Note that this solution shows sin¬ 
gular behavior at rss = 1.804. In the bottom figure (b), the 
conformal structure is depicted. The meaning of the curves 
and symbols in this figure are same as those of FIG. 3(b). 


C. Ultimate thunderbolt singularity (uTS) 

: the case of 32 7 ^ 0 and f5i = P 2 = 0 

Next, we consider only the spatial higher curvature 
correction term Ti? to take into account the back reac¬ 
tion effect of the Lifshitz scaling in the high energy limit. 
That is the case of 32 7 ^ 0 and fix = ^2 = 0. From the 
discussion in section IIIB, the graviton horizon does not 
change and its position is still at rcH where Tc{rQxi) = 0. 
On the other hand, the scalar-graviton horizon is shifted 
from the Killing horizon to the universal horizon because 
of term in the sound speed (3.19). 

We first assume 52 < 0. In Fig. 9(a), one numeri¬ 
cal solution is shown for g 2 = —1. The setting of the 
coupling constants and the boundary conditions are the 
same as those of Fig. 1 except g 2 = —1. Notably, we find 
that a singular behavior on the metric horizon found for 
untuned arbitrary value of 02 in the ae-gravity theory 
vanishes. The metric horizon turns to be regular for any 
value of 02 . Instead, a singular behavior is found inside 
the metric horizon. The {v,r) component of the metric, 
B{r) vanishes there. The aether field aligns perpen¬ 
dicular to the timelike Killing vector near the singular 
point, which corresponds to the universal horizon with 
u ■ ^ = 0. 

We calculate the four-dimensional Kretchmann invari¬ 
ant and the quadratic scalar of the three- 

dimensional Ricci tensor , which are shown in 

Fig. 9(b). Near the singular point, those two scalar 
functions diverge. Thus, this singular point is a phys¬ 
ical singularity rather than a coordinate singularity. So 
we find that the universal horizon becomes singular. Ad¬ 
ditionally, the structure of this spacetime is depicted in 
FIG. 9(c). 

What causes this type of singularity ? In order to 
clarify this question, we first show the relation between 
g 2 and the radii of the graviton horizon rcH, the metric 
horizon tkh and the singular universal horizon run in 
Fig. 10. We also give the detailed data of the singular 
universal horizon radius run near (72 = 0 in Table II. 

From these results, we find that the universal horizon 


fl2/mgi ruH 

0 0.7200 -5.00 X lO""* 

-1.00 X 10-'‘ 0.7783 -6.00 x 10"'‘ 

-2.00 X 10“'‘ 0.7854 -7.00 x lO""* 

-3.00 X 10"'‘ 0.7893 -8.00 x 10"'‘ 

-4.00 X 10“'‘ 0.7922 -9.00 x lO""* 


ruH 

0.7944 

0.7962 

0.7977 

0.7990 

0.8002 


TABLE II: The detailed value of run for near g 2 = 0. The 
coupling constants and the boundary conditions are same as 
Fig. 10 except g 2 - The values of run are evaluated at the 
point where U{r) = —0.025. 


radii change smoothly from 52 = 0 case to (72 7 ^ 0 cases. 

One may wonder whether it can be regular if we tune 
the free parameter 02 just as the scalar-graviton horizon 
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(a) The evolution of each components of metric and aether. 



(b) The corresponding evolution of the curvatures. 


run : thunderbolt i~^ 



(c) The Carter-Penrose diagram 

FIG. 9: The spherically symmetric solution with higher spa¬ 
tial curvature correction. The setting of the coupling con¬ 
stants and the boundary conditions are the same as those of 
Fig. 1 except — —1.000. In the top figure (a), the 

blue, red and green curves indicate the functions T(r), B(r) 
and a{r), respectively. The dashed purple and dashed brown 
curves indicate the Tq and U{r). The graviton horizon is 
found at rcH = 1.248, while the metric horizon exists at 
rKH = 1.250. The calculation has been broken down at 
ruH = 1.157 where U{r) approaches to zero. In the middle 
figure (b), the dotted gray and dotted black curves indicate 
ln(i?^,/po-7?''‘'^'^) and In(77^,/77.'“'), respectively. In the bottom 
figure (c), the conformal structure is depicted. The meaning 
of the curves and symbols in this figure are same as those of 
FIG. 3(b). 


in the ae-gravity theory. To see this, we shall perform the 
expansion of the basic equations around the universal 



FIG. 10: The relation between several horizon radii and 32 - 
The coupling constants and the boundary conditions are the 
same as Figs. 1 and 9 except < 72 ■ The purple, blue and 
brown star marks indicate the graviton horizon radius ron, 
the scalar-graviton horizon radius rsn and the singular uni¬ 
versal horizon radius run in the m-gravity theory, respectively. 
The purple, blue and brown round circles indicate the gravi¬ 
ton horizon radius tgh, the metric horizon radius ckh, and 
the radius of the singular universal horizon run. Numerically 
the position of run is evaluated at the point of U (r) = —0.025, 
because U{r) = 0 is singular. 


horizon. Focusing on the coefficients of the most highest 
r derivative terms, T'"{r), B"'{r) and a"'{r) terms in 
the (0, 9) component of the Einstein equation, we find all 
of them have been vanished at the universal horizon (see 
Appendix.C). This result does not depend on the value of 
02 . This fact means that the singularity on the universal 
horizon cannot be remedied by tuning the free parameter 
02 unlike the ®-case. While there is no singular behavior 
on the scalar-graviton horizon rsn in the infrared limit 
where Ts vanishes. 

Hence it is not quite unnatural to consider that the 
singular behavior which appears on the scalar-graviton 
radius with Ts = 0 in the infrared limit is shifted to the 
universal horizon when we include the higher curvature 
term from the Lifshitz scaling. Namely, the dispersion 
relation of the scalar-graviton with nonzero g 2 gives the 
infinite sound speed of the scalar-graviton. The similar 
situation is found for the exact solution with C14 = 0 in 
the ae-gravity theory [13], for which the sound speed of 
the scalar-graviton (3.12) becomes infinite and the scalar- 
graviton horizon coincides with the universal horizon. In 
this case, however, this singularity can be removed by 
choosing an appropriate value of 02 . 

We shall turn our attention to the physical interpreta¬ 
tion of this solution. Recall that the universal horizon is 
defined by surface where u-^ = 0, i.e., a static limit for the 
ultimately excited dispersive particle with z ^ 1 Lifshitz 
scaling. In other words, only the particle which possesses 
infinite energy can stay on the surface, and any future- 
directed signal on the surface cannot goes outward even if 
the particle is spacelike with infinite energy. Namely, the 
information on the universal horizon must not be leaked 
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out to outside of the horizon. Although this solution can¬ 
not be regarded as a black hole solution whose spacetime 
singularity is isolated by an event horizon, the cosmic 
censorship hypothesis is not violated on this account. 

This singularity is very similar to the thunderbolt sin¬ 
gularity if we replace a null event horizon with a spacelike 
universal horizon. The universal horizon is a real horizon 
for the z > 1 Lifshitz scaling particles. So the singular 
universal horizon cannot be detected by any outside ob¬ 
servers. Although there is no causally disconnected re¬ 
gion, it is not a naked singularity. We then call it a ulti¬ 
mate thunderbolt singularity (uTS). We may speculate 
that the appearance of the thunderbolt singularity on the 
universal horizon via TZ^ term indicates quantum gravi¬ 
tational loop correction in the Lorentz violating system 
as in quantization of the Lorentz invariant system with 
a thunderbolt singularity. 

Besides, one may wonder about solutions for the pos¬ 
itive value of g 2 - In fact, this case is less interest¬ 
ing. Namely, all black hole horizons which exist in the 
case of 52 = 0 completely disappear. The property 
of the solution is quite unphysical, i.e., the functions 
T(r) and B{r) show positive divergence at smaller ra¬ 
dius r ^ 0[\QGnM), while the function a{r) drops to 
zero without forming any horizon. In other words, there 
does not exist any type of black holes discussed before 
(See Fig. II). 



FIG. 11: The typical example of the solution with positive 52 . 
We set 02 = 1.123 X 10~® and 52 /mpi = 1.000. The remaining 
coupling constants and the boundary conditions are set to the 
same values as those of Fig. 9. The blue, red and green curves 
indicate the functions T(r), B(r) and a{r), respectively. The 
dashed purple and dashed brown curves indicate the Tq and 
U{r). The function a{r) drops to zero, while T{r) and B{r) 
diverges for smaller radius r. 


D. Solutions with z = 2 Lifshitz scaling terms 
: the case of 52 / 0 and Pi, 132 / 0 

If we have only term, we find uBH for an appro¬ 
priate value of 02 . The universal horizon is not singular. 
On the other hand, when we have only TZ^ term, the 
universal horizon becomes singular, giving a thunderbolt 


singularity uTS. One may wonder what happens if both 
z = 2 Lifshitz scaling terms, u* and TZ^, exist (52 ^ 0 and 
Pi ^0). Since TZ^ and iPTZ give the highest derivative 
terms in the equation of motion, namely, T"{r), B"(r) 
and a'"{r), the uTS spacetime with \Pi\ <C I 52 I is not 
so different from the original one. On the other hand, 
if 1 , 01 1 is not so small compared with I 52 I, the spacetime 
tends to generate a singularity caused by the aether field. 
More specifically, a singular spherical shell appears be¬ 
fore forming the thunderbolt singularity on the univer¬ 
sal horizon. Thus, we can only find iBH (ii) and tNS 
spacetime in this situation. 

Finally, we mention the case of ,02 0. Although we 

have shown that there must not exist a regular universal 
horizon in this case (see Appendix C), any horizons can¬ 
not be found as far as our numerical analysis. In other 
words, the spacetime produces a singular spherical shell 
before forming any horizon, namely, tNS solution. 


V. PROPERTIES OF SOLUTIONS 

We shall discuss the properties of obtained black hole 
and thunderbolt singularity solutions from several view 
points. 


A. Distribution of the aether field 

In our solutions, there are two free parameters, M 
and 02 . The mass is a conserved quantity, which char¬ 
acterizes the solution. Althogh the different value of 02 
gives the different solution, 02 may not correspond to any 
conserved quantity. In order to understand the physical 
meaning of 02 , we consider the energy density distribu¬ 
tion of the aether field. 

We define the effective energy density and pressures of 
the aether field by 

:= , 
p[r] - 

^ as • ^ ^ ^ 1 

:= + ( 5 . 1 ) 

frl f±l 

When we calculate Pis,PL ,PL we use the Einstein 
equations: 

SttGTI:! = . (5.2) 

[ael 

It makes easy to evaluate once we obtain the solu¬ 
tions. 

H f-Ll 

First we show p^, Pi , and Pi for the ae-black hole 
in Fig. 12 as a reference. In the ae-black hole, all aether 
quantities are always negative. Additionally, the follow¬ 
ing quantity Pga is introduced in order to examine the 
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FIG. 12: The distributions of the aether cloud for the ae- 
black hole. The solid cyan, dashed magenta, dotted yellow 
and dot-dashed black curve indicate p^., and E^., 

respectively. The coupling constants and the boundary con¬ 
ditions are set to the same values as those of FIG. 1. 


strong energy condition, 



It is also negative definite, which means that the strong 
energy condition is broken. 

\r] f_Ll 

We then show p^, Pi , Pi and Eg. for each solution 
in Figs. 13 and 14. 

When we add the z = 2 Lifshitz scaling terms, the 
distribution of the aether field is drastically changed. For 
the case of /3i ^ 0, /32 = 52 = 0 (for iBH, uBH and 

iTS spacetime), the localized aether cloud is formed, i.e., 
frl 

Pg. and Pi are localized near the graviton and the scalar- 
graviton horizons. We may speculate that the positivity 
of the aether density and radial pressure relax a singular 
behavior at the scalar-graviton horizon, which exists in 
as-theory with general value of 02 . Note that the strong 
energy condition is satisfied for some finite radial region. 
Furthermore, referring FIG. 2 and the aether distribution 
FIG. 13 (a), (b) and (c), for the larger value of 02 , the 
more dense aether cloud forms. As a result, the iTS 
spacetime is emerged in the large 02 region due to the 
gravitational collapse of the aether cloud. 

For the case of 32 0, /3i = /?2 = 0 (for uTS space- 

time), the distribution is different. Referring FIG. 14, al¬ 
though the radial pressure is positive, the energy density 
becomes negative. The strong energy condition is broken 
in the whole spacetime. Note that the radial pressure di¬ 
verges where the shell singularity appears. 


B. Preferable Black Holes 

Although black hole thermodynamics in as-gravity the¬ 
ory has been discussed in the last decade, the complete 
understanding have not yet achieved. Hence, in this sec- 





FIG. 13: The distributions of the aether cloud for the solu¬ 
tions which are found in /3i 0, P 2 = 92 = 0. We illustrate 

the typical example of iBH (i) (the top figure (a)), the uBH 
(i) (the middle figure (b)) the iTS (the bottom figure (c)) in¬ 
stead of showing all kind of the solutions. The solid cyan, 
dashed magenta, dotted yellow and dot-dashed black curve 
indicate pg., P^\ Pj.^' and Eg., respectively. The coupling 
constants and the boundary conditions for (a), (b) and (c) 
are set to the same values as those of FIG. 3, 5, and 7, re¬ 
spectively. 


tion, we only discuss which solution is more preferable 
from the view point of the thermodynamical stability 
In the previous section, we find two-parameter “black 
hole” solutions: One free parameter is a black hole mass 
M, which is used to normalize the variables, and the 
other free parameter is 02 . However we have only one 
Noether charge with respect to time translational sym¬ 
metry, which is the black hole mass M given by (3.10). 
As we showed in the previous subsection, the parame¬ 
ter 02 describes the distribution of the aether field, but 
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FIG. 14: The distributions of the aether cloud for the uTS 
spacetime. The solid cyan, dashed magenta, dotted yellow 
and dot-dashed black curve indicate pa., and Pae, 

respectively. The coupling constants and the boundary con¬ 
ditions are set to the same values as those of FIG. 9. 


does not provide a conserved quantity. It just describes 
a cloud of the aether held around a black hole or a thun¬ 
derbolt singularity. 02 describes a different configuration 
of the aether cloud. 

Hence, fixing a black hole mass M and changing 02 , we 
may find most preferable configuration of the aether field, 
which gives a stable solution. To find such a solution, we 
adopt the view point of thermodynamical stability, i.e., 
we assume that the maximum entropy determines the 
stable configuration. 

However, the definition of the black hole entropy is un¬ 
clear due to the unavailability of Wald’s Noether charge 
method on the black hole horizons in ae-theory and its ex¬ 
tension. In addition, according to [38], the black hole en¬ 
tropy is modified by the higher curvature terms. Namely, 
it is given by the integration of functional derivative of 
the action with respect to four dimensional Riemann ten¬ 
sor denoted by over the bifurcation surface B. In 

our case. 


= (1 - /32U^ - 2g2TZ)gPPg''‘^ . (5.4) 

Unfortunately, this modification factor (I — /32U^ — 2p2^) 
diverges near the universal horizon due to the singular 
behavior of the three curvature TZ. 

Hence here we consider only the case without the 
higher-curvature correction terms (32 = P 2 = 0)- Then 
we simply assume the black hole entropy is given by the 
area of the horizon yl(r//), where th is one of horizon 
radii If we find more appropriate definition of the 
black hole entropy, our result would be changed. 


According to [38], the black hole entropy is modified by the 
higher curvature. Namely, it is given by the integration of func¬ 
tional derivative of the action with respect to four dimensional 
Riemann tensor over the bifurcation surface B. When term is 
considered, the entropy should be same as ae-theory’s one. It is 
because never produces any additional term by the functional 
derivative. 


Here we adopt the universal horizon to evaluate the 
black hole entropy: 


dGN Gn 


(5.5) 


In Fig. 15, we show the radii of the universal hori¬ 
zon ruH with respect to 02 for uBH spacetime. The 



FIG. 15: The radii of the universal horizons run with respect 
to 02 for uBH spacetime. The solid red, dotted orange, 
dashed green and dot-dashed blue curve indicate for Pi /rripi = 
0.250, 0.500,1.000 and 2.000 case, respectively. The remaining 
coupling constants and the boundary conditions are set to the 
same values as those of Fig. 1. The points which give the 
largest value of run for each /3i are indicated by the circle 
plot on the each curves. The turning point where the curve 
has the sharp edge is given by a 2 (turn) = 0.373 for every /3i. 


property of the universal horizon is summarized as fol¬ 
lows : Each universal horizon radius is a convex upward 
function with respect to 02 except the right edge. The 
maximum value of black hole entropy for each value of 
Pi is given at the top of the convex, which is denoted 
by ?’uH(max)- At the end of the convex function which 
is called a ’’turning point” denoted by a 2 (turn)i tbe func¬ 
tion sharply bounce and turns to increase until the reg¬ 
ular universal horizon disappears. Beyond this point, we 
find iTS solution. Remarkably, the value of 02 at turn¬ 
ing point and the right edge of these plot are invariant 
with respect to /3i, which are given by a 2 (turn) = 0.373 
and a 2 (end) = 0.377, respectively. Note that a 2 (end) cor¬ 
responds to the border between uBH and iTS space- 
times. 

We give the detailed data of the thermal quantities 
with the maximum entropy and at the turning point for 
each value of /3i in TABLE HI. 

When we take the maximum value of the universal 
horizon radii with respect to 02 , we may find the most 
preferable black hole solution for a given mass M. It is 
because such a maximum point may give a stable solution 
from the view point of the black hole thermodynamics. 

Then we find a series of these most preferable black 
hole solutions in terms of M. In Fig. 16, we plot the 
horizon radii of such a black hole v.s. the gravitational 
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Pi! mil 

®2(max) 

^ UH(max) 

0^2 (turn) 

^UH(turn) 

0.250 

0.296 

0.607 

0.373 

0.345 

0.500 

0.301 

0.499 

0.373 

0.262 

1.000 

0.304 

0.380 

0.373 

0.178 

2.000 

0.307 

0.259 

0.373 

0.103 


TABLE III: the detailed data of the universal horizon radii at 
the maximum point and at the turning point for each value 
of /3i. 


mass M for = 1. Additionally, that of the ae- 

black hole, namely, Pi = P 2 = 92 = ^ case, are also 
shown as reference. For all horizons (the graviton, scalar- 


Pi! mil 

^2(max) 

^UH(max) 

0 

1.112 X 10"“ 

‘ 0.720 

0.001 

0.229 

1.001 

0.005 

0.248 

0.949 

0.010 

0.259 

0.916 

0.050 

0.281 

0.799 

0.100 

0.288 

0.727 


TABLE IV: The detailed values of run and corresponding 
02 of the most preferable uBH solutions for j3i <C 1. The 
case of /9i = 0 corresponds to the ae-black hole. The coupling 
constants and the boundary conditions are set to the same 
values as those of FIG. 1. 



FIG. 16: The mass dependence of the horizon radii for the 
most preferable black hole solutions with /Si/rripj = 1 are in¬ 
dicated by solid curves. 02 is fixed so that the maximum 
universal horizon radii is obtained. The remaining coupling 
constants and the boundary conditions are set to the same 
values as those of Fig. 1. The brown, purple and blue curves 
indicate the radii of the universal horizon, the graviton hori¬ 
zon and the scalar-graviton horizon, respectively. The radius 
of the uBH universal horizon turn to be greater than that 
of m-black hole in M > Merit = 1.497mpi region. 

graviton and universal horizons) of the ae-black hole, their 
radii seem to increase linearly in terms of the mass M, 
which is the same as that of Schwarzschild solution. 

Whereas, for our uBH, the universal horizon radius 
increases with a higher power-law function of M than the 
linear one in the small mass region, and it approaches a 
linear function for large values of M. The uBH univer¬ 
sal horizon radius is smaller than ®-black holes’s one for 
M > Merit = 1.497TOpi. We shall refer Merit as a critical 
mass. Note that the universal horizon radius seems to 
vanish for small value of M, but it is not clear whether 
it vanishes at a finite mass Mmin(> 0), or at M = 0. 

One may wonder whether our uBH solution will re¬ 
cover the ae-black hole or not when /3i approaches to zero. 
To see this, we give the detailed data of ruH(max) and cor¬ 
responding a 2 (niax) near /3i = 0 in Table IV. From this 
table, we find that the value of gradually decreases 

as /3i —>■ 0. Whereas, the maximum value of the univer¬ 


sal horizon radius increases in such a limit. Thus, we 
can conclude that the uBH solution cannot be smoothly 
connected to the as-black hole. 


C. Smarr’s formula and Black Hole Temperature 


Concerning the black hole first law in ^-theory, it 
is found that the aether held prevents from estab¬ 
lishing the black hole mass-entropy relation on the 
Killing horizon[31] via Noether charge method[38, 39]. 
Notwithstanding, the Smarr’s formula in ^-theory has 
been proposed only in static and spherically symmet¬ 
ric conhguration[13], which is established by applying 
Gauss’s law to the aether held equation. It is found that 
the aether portion of the basic equation can be reduced 
to Maxwell-like form in static and spherically symmetric 
spacetime : 

= (5.6) 


where, is given by 


Ffiy :— , 

■Ci4 


(5.7) 


q := 


1 2 


- Ci3 


(u • s)(u ■ C) + (1 - C13 )k 


+ £|i(V.u)(s.^), 


(5.8) 


s^ is a spacelike unit vector perpendicular to and a 
surface gravity k, which is given by is 

equivalent to the following form due to a spacetime sym¬ 
metry : 


K = 



(V„^^)(V“e^). 


(5.9) 


Note that the structure of (5.6) is similar to Maxwell’s 
equation, thus, we can perform the hux integration. Ac¬ 
cording to Gauss’s law, the integration of over Br 
which is a two sphere at radius r must produce same 
value for any r. Considering the integration over spatial 
inhnity Boo and Br, then we obtain 


M = q{r) 


47rGN ’ 


(5.10) 
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where A{r) := AttA is a surface area of Br- Note that 
this relation is held for any r. When we evaluate the 
r.h.s. of Eq. (5.10), we find the Smarr’s formula: 

M = , (5.11) 

where = ^(^uh)/'?’' is a black hole temperature. 

Turning to our attention to the case of including higher 
curvature and aether effects, although the Maxwell-like 
aether equation is not found, we may obtain the relation 
between the black hole mass and the entropy (or the area 
of the universal horizon) by evaluating the deviation from 
the ae-black hole’s one. Note that our thermodynamical 
analysis is limited to the black hole solution only with /3i 
term. If we consider the surface gravity or black hole tem¬ 
perature on the singular horizon, these quantities must 
diverge. Thus, we shall not discuss the solutions without 
a regular universal horizon. 

We shall presume the mass-entropy relation as follows: 

M (5.12) 

where := T^bh + introduced 

as the correction by the z = 2 Lifshitz scaling term from 
the ae-black hole temperature. 

In Fig. 17, we show the mass-temperature relation of 
the most preferable black hole solutions defined in the 
previous subsection. We also show 8e-case[13] as a refer¬ 
ence. From this plot, we find that the ®-black hole tem- 



FIG. 17: The mass-temperature on the universal horizon for 
the most preferable black hole solntions. The solid red, dot¬ 
ted orange, dashed green and dot-dashed blue curves indicate 
the temperatures foi' /3i/^pi = 0.250, 0.500, 1.000 and 

2.00, respectively. The values of 02 for each case are fixed so 
that the maximum universal horizon radii is obtained. The 
remaining coupling constants and the boundary conditions 
are set to the same values as those of Fig. 1. The m-case is 
also shown as a reference, which is indicated by the two-dot 
chain black line. The critical mass for each j3i are indicated 
by the colored dot. 

perature seems to be inversely proportional to the black 
hole mass M, which is the same as the Schwarzschild 


black hole in GR. Whereas, that of our uBH solution 
is clearly far from the inverse M law. More specifically, 
the black hole temperature of uBH turns to be lower 
than that of ae-case for the range of M > Merit, where 
Merit is the critical mass defined in the previous subsec¬ 
tion. In TABLE. V, we show the values of Merit and 
corresponding ryn in terms of /3i. 


/3i/mpi Merit /m-pi run/mpi 


1.000 

1.497 

2.134 

0.900 

1.420 

2.029 

0.800 

1.339 

1.913 

0.700 

1.253 

1.789 

0.600 

1.160 

1.656 

0.500 

1.059 

1.512 

0.400 

0.947 

1.352 

0.300 

0.820 

1.171 

0.200 

0.670 

0.956 

0.100 

0.473 

0.676 


TABLE V: The critical masses and corresponding universal 
horizon radii in terms of /3i. 


From FIG. 17, the temperatures of the uBH solutions 
seem to obey the inverse M law in large M region, while 
it may decrease exponentially in small M region. To 
examine this behavior, we shall assume that the mass- 
temperature relation is given by the following functional 
form : 

T^bh'^ « M - 1) 

f for M — Mniin <C v 

y^iv/ (M - Mjnin) for M - Mjnin » V ’ 

(5.13) 

where n, v and Mmin are some fitting parameters. Since 
the universal horizon radius is related to the black hole 
mass as (5.12), the universal horizon should disappear 

when M = Mmin, where T^bh diverges. 

The curves in FIG. 17 are approximately reproduced 
by the above functional form (5.13) if we choose the fit¬ 
ting parameters given in TABLE. VI. In our numerical 
result, it is found that Mmin is extremely small, or pos¬ 
sibly vanishes. Moreover, it is notable that the values of 
[IV seem to be universal for j3i > 0, namely, all asymp¬ 
totic behaviors in large M region, which is denoted by 
^ TOpi), are the same, but differs from that of 
the ae-black hole. 


VI. CONCLUSION AND DISCUSSION 

Without Lorentz symmetry, the Killing horizon is no 
longer the event horizon in the static and spherically 
symmetric spacetime due to the presence of the super¬ 
luminal propagating modes. However, in the context of 
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dl/m-pl 

^ TTlpl 

I'lmpi 



0.250 

0.136 

0.649 

6.561 X 10"' 

0.088/M 

0.500 

0.096 

0.919 

9.261 X 10"^ 

0.088/M 

1.000 

0.068 

1.299 

9.468 X 10"^ 

0.088/M 

2.000 

0.048 

1.837 

9.998 X 10"^ 

0.088/M 

0 

N/A 

N/A 

N/A 

= 0.151/M 


TABLE VI: The fitting parameters in asymp¬ 

totic behaviors in large M region for each value of j3i. The 
coupling constants and the boundary conditions are set to the 
same values as those of FIG. 17. We also show the temper¬ 
ature of the ae-black hole as a reference in the bottom 
line. 


the ^-theory or infrared limit of the non-projectable HL 
gravity, a black hole solution can be still constructed by 
considering the universal horizon which is a static limit 
for an instantaneously propagating particle. In this pa¬ 
per, we have studied the backreaction to the black hole 
solution in ^-theory by the Lifshitz scaling terms. We 
have analyzed the ultraviolet modification of the ®-black 
holes including the simple scalar terms with z = 2 Lif¬ 
shitz scaling, specifically, the quadratic term of spacial 
curvature associated with the hypersurface orthogonal 
aether field (7?.^), the quartic term of the aether accel¬ 
eration and the product of the spacial curvature 

and the quadratic term of the aether acceleration {u^TZ). 

Only for the case with term (/3i ^ 0 and /32 = 
gi = 0 ), we have succeeded in finding a black hole so¬ 
lution with regular universal horizon, which is referred 
as uBH. In contrast to the as-case, the black hole so¬ 
lutions are obtained without tuning the boundary pa¬ 
rameter of the aether field 02 . In other words, uBH 
solutions exist in a finite range in the (/3i, 02 ) parameter 
plane (see FIG. 2). If we select the parameter 02 be¬ 
yond this range, a singular shell appears before forming 
a universal horizon. While, considering u^TZ and/or TZ^ 
terms (/32 7 ^ 0 and/or g 2 7 ^ 0 ), any black hole solution 
with regular universal horizon cannot be constructed due 
to the divergence of the basic equations on the universal 
horizon (see Appendix.C). However, including only TZ^ 
term with negative g 2 , we have found the solution with a 
thunderbolt singularity {uTS), whose universal horizon 
turns is still singular but the singularity is not observed 
any outside observers. Although this solution cannot be 
regarded as a black hole, the cosmic censorship hypoth¬ 
esis is not violated. Since the thunderbolt singularity 
had been discovered in the context of the quantum grav¬ 
ity in lower dimension, the emergence of this solution is 
not so strange since the z = 2 Lifshitz scaling terms are 
regarded as quantum gravitational corrections. 

We then have studied several properties of our solu¬ 
tions including spacetime structures with their Carter- 
Penrose diagrams. To investigate the physical meaning 
of 02 , we have shown the effective energy momentum 

[asl 

tensor of the aether field , i.e., the effective energy 


density and pressure. From these results, the parameter 
02 seems closely related to the distribution of the spher¬ 
ical aether cloud. Namely, large value of 02 configures a 
dense aether cloud and eventually induces a collapse of 
the aether field, which results in a formation of a ’’time¬ 
like” singular shell (tNS). Moreover, we may speculate 
that regularity on the (infrared limit of) scalar graviton 
horizon is recovered due to the localization of the aether 
field near the horizon, which does not appeared in ae-case. 

Finally, we have explored the thermodynamical aspect 
of the uBH solution. Since the maximum universal hori¬ 
zon radius is obtained by choosing an appropriate value 
of 02 for a given mass M, the parameter 02 may be fixed 
so that the area of the black hole, which may be regarded 
as the black hole entropy, becomes maximum. This so¬ 
lution may provide the most preferable black hole be¬ 
cause it can be thermodynamical stable. Additionally, 
the Smarr’s formula and the black hole temperature are 
also examined. When we speculate the mass-temperature 
relation as (5.12), it is found that the temperature does 
not obey the inverse M law at least for small M unlike 
ae-case. It increases exponentially as M approaches to 
Almin, for which the horizon area vanishes. 

Although there are several coupling constants in the 
action we consider (2.5), we have manipulated only ul¬ 
traviolet coupling constants, i.e., /3i, P 2 and g 2 from z = 2 
Lifshitz scaling terms. The coupling constant C 13 , C 2 and 
Ci 4 are all fixed so that the ae-black hole is restored in 
infrared limit. However, we shall emphasize that the ex¬ 
istence of the other types of black hole solution should 
not be excluded in this theory with the different values 
of the coupling constants C1-C4. For example, non-black 
hole solution in ae-theory (a static and spherically sym¬ 
metric solution with Killing horizon but without a uni¬ 
versal horizon found in [30]) may turn to form a causal 
boundary due to the Lifshitz scaling terms. 

Turning our attention to more energetic region, it is 
obvious that the z = 3 Lifshitz scaling terms which are 
required by power-counting renormalizability of gravity 
turn to be dominant rather than z = 2 terms. Namely, 
our analysis in this paper should correspond to the in¬ 
termediate region between the energy scale described in 
ae-theory and the such a ultimately high energy scale. If 
we consider more ultraviolet modification terms in ac¬ 
tion, we may obtain the additional intriguing spacetime 
such as a singularity-free solution which is discovered in 
the context the early universe in HL theory. 


ACKNOWLEDGEMENT 

YM would like to thank T. Wiseman for giving an op¬ 
portunity to start this study and useful discussion when 
he was staying at Imperial College London as the Eras¬ 
mus Mundus PhD fellow. He also thanks to T. Kitamura 
for useful comments. KM acknowledges S. Mukohyama, 
N. Ohta, and S.M. Sibiryakov for discussions in the early 
stage of the present study. This work is supported by a 









20 


Waseda University Grant for Special Research Projects 
(project number : 2015S-083) and by Grants-in-Aid from 
the Scientific Research Fund of the Japan Society for the 
Promotion of Science (No. 25400276). 


Appendix A: The disformal transformation in a 
static and spherically symmetric spacetime with 
asymptotic flatness. 

In this section, we shall show the transformation law in 
Eddington-Finkelstein like ansatz (3.1) under the disfor¬ 
mal transformation. We would also like to confirm that 
the spacetime properties, i.e., time independence , spher¬ 
ically symmetry and asymptotically flatness are held. 
The disformal transformation we consider is 


9iiu = + (1 - (^)u^Uv , = a . (Al) 


where the original metric and aether are given by 
(3.1). Then, the each components of and is given 
by 


ds^ = —T{r)dv'^ + 2B{r)dvdr + f{r)dr'^ + r'^dfl^ , 

(A2) 

= (a{r),b{r),0,0^ , (A3) 

where. 


f(r) 

B{r) 

f{r) 

a(r) 

S(r) 


r(r) - (1 - cr) 

1 


1 -I- a{r)'^T{r) 


2a(r) 


(A4) 


B{r) [l + a — {1 — a)a{r)‘^T{r)] , (A5) 


(1 — a)a{r)^B{r)^ 


a ^'-a{r), 


a{r)^T{r) — 1 
2a^/^a{r)B(r) 


(A6) 

(A7) 

(A8) 


From these form, we find that the time independence and 
spherical symmetry are held after transformation. Hence 
the horizon radius for the propagating degree of freedom 
whose sound speed is unity in this frame is given by a 
null surface T{r) = 0. 

Note that the (r, r) component of the metric grr = f{r) 
is generated unlike original Eddington-Finkelstein type 
ansatz. It is not quite unnatural considering the geo¬ 
metrical meaning of the disformal transformation. More 
specifically, the transformation can be interpreted as a 
rescaling of timelike separation between two spacelike 
hypersurface with fixing three-dimensional space. Then, 
the light cone whose opening angle is 90° in original {g, u) 
frame is distorted after disformal transformation. Hence, 
the null coordinate v in original [g, u) frame is no longer 
null in ( 5 ,m) frame. 

We then introduce new coordinate system {v *, r) in 
which V* becomes a null coordinate, v* is defined by 


* = 


= 


B{r)-^B{rf+T{r)f{r) 
dv -^- dr 


dv • 


T{r) 

2(1 — B{r) 

1-b fTi/2 _ (1 _ CTi/2)a(r)2T(r) 


dr 


(A9) 


dv 
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Then, the metric and the aether field are transformed 
into 


df = -f*{r)dv*'^ + 2B*{r)dv*dr + r'^dVt^ , 

= (^d*{r),b*{r),0,0^ , (AlO) 

where, 


f*{r) 

B*{r) 

d*{r) 

b*{r) 


a ^T{r) + (l - CT 


1 + a{r)'^T{r) 
2a{r) 


B{r), 


(All) 

(A12) 


2a^/^a{r) 

1 + crl/2 _ (1 _ crl/2)a(,.)2r(r) ’ 

a{ryT{r) — 1 

2a^/'^a{r)B{r) 


(A13) 

(A14) 


Since Eddington-Finkelstein type metric is restored, i.e., 
the (r, r) component of the metric is vanished, we can 
regard v* as a null coordinate in {g,u) frame. 

We now focus on the asymptotic property. Substitut¬ 
ing (3.8), we find 


f*{r) 

= 1 + —+ 0 (r- 2 ), 
r 

(A15) 

B*{r) 

= l + 0 (r- 2 ). 

(A16) 

d*{r) 


(A17) 

b*{r) 

= 0 (r- 2 ). 

(A18) 


Thus, it is found that the asymptotic flatness is held 
even if the disformal transformation is performed. Ad¬ 
ditionally, the mass of the spherical object, i.e., Noether 
charge with respect to time translational symmetry is 
also invariant under the disformal transformation. 


Appendix B: The invariance of spatial curvature 
under disformal transformation 


As we mentioned, the action of ae-theory has an in¬ 
variance except each coupling constants. In this section, 
we shall show the transformation law in detail. For con¬ 
venience, we define a new tensoral quantity as a 

change of Christoffel symbol, 


r Ct _ -pQ! I -yCX. 

M ~ ^ M ^ ^ M ■ 


X' 


Pi 


:= (1-cr) 


— 1 O/. • Ol 

G U /C/3<y — U 


■ (Bl) 


where, := 7 (^“VaUi/) is an extrinsic curvature asso¬ 
ciated with the aether. Note that is imposed 

by the hypersurface orthogonality of the aether. Then, 
covariant derivative of the aether is given by 

VaU/3 = \yaUp - (l - Cr“^) ICafi] , (B2) 


and, we find the aether is invariant, 

= it X . (B3) 

thus, the Ci 3 , C 2 and ci 4 terms in the action (2.5) are 
transformed into 

(V„m'5)(V;3u“) = a-i(V„n'5)(V^u“), 

(V-u)2 = (J-l(V•^i)^ 

. (B4) 

The four dimensional Riemann tensor is trans¬ 

formed into 


= Kp^ + V/3A% - V.A“ ^ 

, (B5) 

and thus, we find the transformation law of Ricci scalar 
R as follows : 

R = i?_(l_cr-l) [{XaU^){XfiU°‘)-{V-uf] 

-I-(total derivative term), (B 6 ) 


Since the total derivative term in the action can be inte¬ 
grated out, we shall abbreviate it. 

Turning our attention to the transformation law of the 
spatial three-curvature From the Gauss-Codazzi 

relation, the spatial curvature can be expressed in terms 
of 5 ^ 1 / and : 

, (B7) 

Since the three metric 7 ^;^ is invariant, we have only to 
consider the terms associated with in (B5) to see 
the transformation of the first term. Then, 

= - (1 - a-i), 

= - (1 - a-i), 

7V.^7"a^“,^V = 0, 

= 0, (B 8 ) 

The last two terms of RHS in (B7) are transformed into 

i^OL /T~^1C IC^ 

^flOC^ y G y ? 

(B9) 


Thus, we can find TZ^i, is invariant under the transfor¬ 
mation and obtain (2.15). 

Since TZfj,i, and are invariant under the disformal 
transformation which is a rescaling of timelike separation 
between two spacelike hypersurface, we expect that its 
spatial covariant derivatives are also invariant. In order 
to confirm it, we examine the transformation law of a 
purely spatial tensoral quantity so that = 

Z^y... and u'^Zau- - = u°‘Z^a. = • • • = 0. Then, 


'DfJ.Ziyp.. 


— d ^d ^ ' 

I fj, lu Ip 

— ty ^y . 

I fJ, lu Ip 


= V^Zyp. 


■ Vr.Z 


a-^/37--- 

XaZp^... -I- X^^pZs^ 


X^^Zps- 


(BIO) 
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The terms proportional to in the second line are 

vanished due to the spatial property of Thus, we 

conclude that all of the spatial higher derivative terms 
derived by non-projectable HL gravity are invariant un¬ 
der disformal transformation. 


Appendix C: The regularity conditions of horizon 

We illustrate the detail of the regularity on the black 
hole horizons with/without z ^ 1 Lifshitz scaling terms. 
In order to investigate the behavior near horizon, we fo¬ 
cus on the highest r-derivative terms in the basic equa¬ 
tions. 


1 . Einstein-aether case : Pi = P2 = 02 = 0 

In this case, the set of the evolution equations, {v,v) 
and {6,0) components of (2.29) and component of 
(2.30), can be simplified into following form : 



= T''[T,T',B,a,a'], 

(Cl) 

B' 

= B'[T,T',B,a,a'], 

(C2) 

a^' 

= a"[T,T',B,a,a'], 

(C3) 


To see the cause of singular behavior on the scalar- 
graviton horizon, we expand (C1)-(C3) around rsn = 0. 
We then find 

T"{r) = M^ + ^T[„](rsH)[Ts(r)]-, (C4) 
Ts{r) 

B'ir) = ^^^^ + i2^H{rsH)[Ts{r)r, (C5) 

a'{r) = ^^^^T7^ + Xl“N(^SH)['rs(0]”- (C6) 

I 


where rs(r) = 7s(rsH)(?’-^ sh) and T[^], a[a.], T[„], 

i3[„] and a[„] are functionals with respect to r'(rsH), 
S(?’Sh), n('rsH) and a'(rsH)- Obviously, the irregularity 
on the scalar-graviton horizon is due to the coefficient of 
Ts{r)~^ terms. As a result, for regularity, all of 
and a[ae] must vanish, otherwise the evolution equations 
diverge at r = rsn- Although the explicit form of the 
coefficients of Ts{r)~^ terms are quite complicated for 
general setting of the coupling constants C 13 , C 2 and C 14 , 
it can be reduced to slightly simplified form considering 
the disformal transformation (2.12). By the disformal 
transformation with a = Cg, we can set the sound speed 
of the scalar graviton to be unity without loss of gener¬ 
ality. In this frame, the scalar-graviton horizon coincides 
with the metric horizon, i.e., Ts(rsH) = T(?'sh) = 0. 

Transforming into the Cg = 1 frame reduces the 
three-dimensional parameter space of coupling constants 
(ci 3 ,C 2 ,ci 4 ) into two-dimensional one. In other words, 
one of these coupling constants is expressed by the other 
two. More explicitly, for example, we can eliminate C 2 
from the evolution equations as follows: Performing sim¬ 
ple calculation, we find Cg = 1 for any C 13 and C 14 , if we 
set 


—2ci3 + 2ci4 — 

2 — 4ci4 -I- 3ci3Ci4 


In this frame, we finally obtain the explicit forms of Tiaj], 
and a[a.] as 


^[ae] 

Bias] 


' 2ci3 - 2ci4 + C13C14 - 2(1 - Ci3)ci4T'(rsH)g^(?^SH)>’SH 
2 ci 4(1 — ci 3 )i?(rsH)a^(rsH)rsH 

_Cl4-S(rSH)_ 

8r(ci4 - 2)(2 - 4ci4 + 3ci3Ci4)a4(rsH) 


-S[33_„]a”(rsH) 


Bi 




(2 — ci 3)(2 — ci 4 )a(rsH) — 2(1 — 013)0140'(r'sH)?^SH 

2ci 4(1 — Ci3)i?(rsH)rsH 


B!, 


(C8) 

(C9) 

(CIO) 
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where 


■B[»,o] := [4 ci 3 - 4(1+ ci3 - c?3)ci 4 + (4-3ci3)cy a'(rsH)?'sH , 

'■= 8(1 - ci3)^ci4a'(rsH)?'SH , 

-^[33,2] := — 4(1 — Ci3)(2 — 2ci4 + C13C14) , 

^[as,3] ■= 2 [—8 + 4 ci 3 + 4(2 — Ci3)^Ci4 — (4 — 3 ci3)c 44] T'(rgH)a^(7’SH)''’SH > 

B[!e,A] ■= “ Ci3)^Ci4T'(rsH)?'SH , 

B[2e,5] ■“ 0 I 

-B[aB,6] := [4 ci3 - 4(1+ C13 - Ci3)ci4 + (4-3ci3)ci 4] r'^(rsH)rsH , (Cll) 


Thus, we conclude that = 0 should be imposed for 
regularity of the scalar-graviton horizon. 

As for the regularity on the universal horizon, we 
also expand the basic equations for T, B and a around 
the universal horizon run, where (it • ^)(r'uH) = 0. 
We then find there is no term with negative power 
of U{r) := U'{r\jH){r — run)- This means that the 
universal horizon is always regular, if it exists, unlike 
the scalar-graviton horizon. 


the scalar-graviton horizon, it is always regular without 
tuning 02 . Similarly, it can be confirmed that the univer¬ 
sal horizon is always regular in the same way. Namely, 
the black hole solutions turn to depend on the mass pa¬ 
rameter Ti as well as the extra parameter 02 . Actually 
we find the black hole solutions in a certain range of 02 
in Section IV B. 

3. The case with TZ^ and/or ii^TZ term : 

(/2 / 0 and/or 132 ^ 0 


2. The case with term : /3i ^ 0 and g 2 = ^2 = 0 

The set of the evolution equations can be decom¬ 
posed into T{r), B(r) and a(r) equations similar to the 
Einstein-aether’s case (C1)-(C3), i.e., the linear-order dif¬ 
ferential equation for B{r) and the second-order differen¬ 
tial equations for T{r) and a{r). The expanded equations 
around rgn = 0 are given by 

r"(r) = ^T(„^^)(rsH)[Ts(r)]", (C12) 

n>0 

B'{r) = ^Bfff^(rsH)[Ts(rr, (C13) 

n>0 

= E<]^(^sH)[rs(r)]r (C14) 

n>0 

where Ts(r) = rg(rsH)(?' - J'sh) and and 

are functionals with respect to r'(rsH), B(rsH), 
^(^Sh) and a'(rsH)- Since all of these equations have no 
terms with negative power of Tg (^) in the expansion near 


In this case we find that T"'{r), B"'{r) and a"'{r), 
which are the highest r-derivatives in the equations, ap¬ 
pear only in the (9,0) component of (2.29). This means 
that the evolution equations cannot be separated into 
T(r), B{r) and a{r) equations unlike the Einstein-aether 
only with the term. Therefore, we just focus on the 
coefficients of T"'{r), B"'{r) and a"'{r) terms in the 
{9,9) component of (2.29). To see the behavior near 
the universal horizon where u ■ ^ = 0, we shall express 
these terms using U{r) := U'{r\jii){r — run) instead of 
T(r) := T'(rgH)(r - rgn): 

0 0TT'"(r) + QBB'"{r) + Qaa"'{r) , (C15) 

where 0 t, 0s and 0a are the functional which are given 
by 

^ _ a(ruH)U^(r)l2g2U(r) - ^2ruHU'(ruH)] 

r^uBHruH) 

^ 4g2UHr) 

^ • r^nB^run) ’ 

and 


2U^(r)ll + a(ruH)U(r)][2g2l/(r) - /32WHl7'(ruH)] 
ruH^^(^UH)B^(ruH) 


Clearly, it is impossible to avoid the divergence of this ficients 0 t, 0 s and 0 a must not vanish simultaneously 
equation on the universal horizon because all of the coef- 
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at ruH for regularity. Therefore, we conclude that the 
universal horizon is always singular. As a result, the 


thunderbolt singularity appears if H? and/or ii^TZ terms 
are joined. 
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